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We present a new method to calculate optical properties of strongly correlated systems. It is based 
on dynamical mean-field theory and it uses as an input realistic electronic structure obtained by 
local density functional calculations. Numerically tractable equations for optical conductivity, which 
show a correct non-interacting limit, are derived. Illustration of the method is given by computing 
optical properties of the doped Mott insulator Lai-^Sr^TiOa. 



I. INTRODUCTION 

Optical spectral functions such as conductivity or re- 
flectivity are very important characteristics of solids 
which give us a direct probe of their electronic struc- 
ture. In the past, very powerful numerical techniques p] 
based on density functional theory (DFT) and local den- 
sity approximation (LDA) have been developed, which 
allowed to access the one-electron spectrum in real ma- 
terials via association of LDA energy bands with the real 
excitation energies. This approach works well for weakly 
correlated systems, where, for example, optical proper- 
ties can be directly computed via the knowledge of 
the band structure and the dipole matrix elements of the 
material. Furthermore, for weakly correlated materials 
LDA is a good starting point for adding perturbative cor- 
rections in the screened Coulomb interactions following 
the GW approach Q. 

Unfortunately, the treatment of materials with strong 
electronic correlations is not possible within this frame- 
work. Strong on-site Coulomb repulsion modifies the 
one-electron spectrum via appearance of satellites, Hub- 
bard bands, strongly renormalized Kondo-like states, 
etc., which are no longer obtainable using static mean- 
field theories such as Hartree-Fock theory or LDA. The 
wave functions in strongly correlated systems are not rep- 
resentable by single- Slater determinants and dynamical 
self-energy effects become important, thus requiring a 
new theoretical treatment based on the dynamical mean- 
field theory (DMFT) Q . Recent advances [f| in merging 
the DMFT with realistic LDA based electronic structure 
calculations have already led to solving such long stand- 
ing problems as, e.g., temperature dependent magnetism 
of Fe and Ni Q, volume collapse in Ce 0], and huge 
volume expansion of Pu § . 

In the present work we develop a new approach 
which allows us to calculate the optical properties of 
strongly correlated materials within the combined LDA 
and DMFT framework. We discuss the expressions for 
optical conductivity using self-energies and local Green's 
functions, which are numerically tractable and correctly 
reproduce the limit of non-interacting electrons. We also 
check the limit of strong correlations by applying the 
method to three-band Hubbard Hamiltonian. Results of 



this test reproduce the available experimental and the- 
oretical data with very good accuracy. We demonstrate 
the applicability of the present scheme on the example 
of doped Mott insulator Lai_ x Sr x Ti03, where we com- 
pare the results of our new calculations with the LDA 
predictions and experiment. 

The paper is organized as follows. In the next section 
m we describe the method for calculation of the opti- 
cal conductivity. Application of the method to doped 
Lai-^Sr^TiOs is described and analyzed in section IIIII 
which is followed by conclusions presented in section ITvl 
Some technical details of the calculations and the down- 
folding and upfolding procedures are given in Appendix 

El 



II. METHOD 

To calculate the optical response functions we utilize 
the dynamical mean field approach where the self-energy 
of the many-body problem is approximated by a local op- 
erator which is, however, frequency dependent. A 
physical transparent description of this method can be 
achieved by introducing an interacting analog of Kohn- 
Sham particles, ip^fr, to) = ipkjun which reproduce the 
local portion of the Green's function in a similar way 
as the non-interacting Kohn-Sham particles ^iy(r) re- 
produce the density of the solid in its ground state. This 
spectral density functional approach 9] has an advantage 
that the k-integrated excitation properties (such, e.g., 
densities of states) can now be associated with the real 
one-electron spectra. The optical transitions between the 
interacting quasiparticles V'kju allow the excitations be- 
tween incoherent and coherent parts of the spectra (e.g., 
between Hubbard and quasiparticle bands) which are in- 
trinsically missing in static mean-field approaches such 
as DFT but are present in real strongly correlated situa- 
tions. 

In order to find the quasiparticles living at a given 
frequency u> we solve the Dyson equation with the LDA 
potential V e ff and the frequency dependent correction 
T,(u>) - E dc , i.e. 

(-V 2 + V eff + E(w) - S rfc - e k ^X, w - 0. (1) 
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A double counting term E^c appears here to account for 
the fact that V e f / is the average field which acts on both 
heavy (localized) and light (itinerant) electrons. Note, 
that due to non-Hermitian nature of the problem, both 
"right" i/i^and "left" ip L eigenvectors should be consid- 
ered, the latter being the solution of the same Dyson 
equation Q with ip placed on the left. The local Green's 
function is constructed from the eigenvectors and eigen- 
values in the following way 



G{w) 



W + fJ, - €kj u 



(2) 



The local self-energy is calculated from the correspond- 
ing impurity problem which is defined by the DMFT self- 
consistency condition 

G(u>) = (w- E imp - £(cj) - A imp (w)) _1 , (3) 

where Ai mp is the impurity hybridization matrix and 
Eimp are the impurity levels. From known Ai mp (u), 
Ei m p and Coulomb interaction U, the solution of the 
Anderson impurity problem then delivers the local self 
energy S(w). The system of equations {Q, J2J and 
(|3*|). together with an impurity solver, i.e., a functional 
S[A imp (o;), E imp , U], is thus closed. 

Solution of the Anderson impurity model can be car- 
ried out by available many-body technique M such as 
the Quantum Monte Carlo (QMC) method [l(| which 
will be used in our work. In practice @, we utilize 
the LDA+DMFT approximation and treat only the d— 
electrons of Ti as strongly correlated thus requiring full 
energy resolution. All other electrons are assumed to be 
well described by the LDA. The Dyson equation is solved 
on the Matsubara axis for a finite set of imaginary fre- 
quencies iuj n using a localized orbital representation such, 
e.g., as linear muffin-tin orbitals (LMTOs) for the 
eigenvectors ipkjuj- 

The optical conductivity can be expressed via equilib- 
rium state current-current correlation function [12 ] and 
is given by: 

<Tp V (w) = ire 2 / ds(j) l j, J/ (e + uj/2,e- w/2) x 



f{e-w/2)-f{e + <j/2) 



(4) 



where e is free electron charge, f(e) is the Fermi function 
and the transport function 4> p ^(e,e') is defined as 



with V being the unit cell volume and 



(6) 



is expressed via retarded one-particle Green's function, 
Gkj (e) ■ Using the solutions Ckj^ and ipkj^ of the Dyson 
equation Q we express the optical conductivity in the 
form: 



cr^(cj) = - 



4-7T 



g'=±l kjj'. 



de- 



kjj' 



kje kj's+ 

f(e-)-f(e+) 



[1-6 



(7) 



where we have denoted = £±u/2, and used the short- 



cut notations ejK = e kje , e kje = ■ 



kje' 



The matrix elements M k jj' are generalizations of the 
standard dipole allowed transition probabilities which are 
now defined with the right and left solutions i/; fl and ip L 
of the Dyson equation: 

M™' j r(e,e')= (8) 



ip?^ and as- 



where we denoted ip£ jE = ip£ je , i/j kje - Vkje 

sumed that W& e ) + = and (V^ s )~ = V£* £ - Ex- 

pressions (J7J), (0 represent generalization of the optical 
conductivity formula for the case of strongly correlated 
systems, and involve the extra internal frequency integral 
appearing in Eq. J7]). 

Let us consider the non-interacting limit when S(cj) — 
Srf c — * i"f — ► 0. In this case, the eigenvalues e k je = 
tkj+iliiPkjs = l k J>- ^tje = l k ^)* = ( k il and tne matrix 
elements '(e, e') are all expressed via the standard 

dipole transitions |(kj|V|kj')| 2 . Working out the energy 
denominators in the expression (JJJ) in the limit ij — > 
and for leads us to the usual form for the conduc- 

tivity which for its interband contribution can be written 
as: 

<7 M „M = — V (kj\v^f)(kj'\v v \kj) x 



kj'#J 

[/(e k i) - f(ekj>)]6(e kj - e kj , 



(9) 



To evaluate the expression ay,,, (a;) in Eq. JJJ) numer- 
ically, we need to perform integration over e and pay 
a special attention to the energy denominator l/{ui + 
e kjE - — e kj'£+)- ^° calculate the integral over e we di- 
vide frequency domain into discrete set of points £j and 
assume that the eigenvalues e k j e and eigenvectors ipk.j e 
to zeroth order can be approximated by their values at 
the middle between each pair of points. In this way, 
the integral is replaced by the discrete sum over internal 
grid Ei defined for each frequency oj. To deal with the 
strong momentum dependence of l/(u> + £ k je- ~ e k.j'£+)' 
linearization of the denominator with respect to k should 
be performed as it is done in the tetrahedron method of 
Lambin and Vigneron O n the other hand, the differ- 
ence between single poles (expression in square brackets 
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of Eq.lJTJ)), after integration over frequency becomes a 
smooth function of k and can be treated together with 
the current matrix elements, i.e by linearizing the nu- 
merator. The described procedure produces a fast and 
accurate algorithm for evaluating the optical response 
functions of a strongly correlated material. 



III. APPLICATION OF THE METHOD 

To illustrate the method of the optical conductivity 
calculation in a strongly correlated system we chose para- 
magnetic doped Mott insulator Lai-^Sr^TiOa. LDA 
cannot reproduce insulating behavior of this system al- 
ready at x = 0, which emphasizes the importance of cor- 
relation effects. Upon doping the system becomes a cor- 
related metal, which at x = 1 (SrTiOa) should be consid- 
ered as a standard band insulator. Photoemission exper- 
iments ^3 as a function of doping display both a lower 
Hubbard band located at near energies 2 eV below the 
Fermi level Ep and a quasiparticle band centered at Ep. 
Previous DMFT based calculations HE! of the density 
of states used t 2g degenerate bands of Ti found near Ep 
and reproduced both these features with a good accuracy. 
The studies of the optical properties for LaTiC>3 with the 
less accurate LDA+U method have been also carried 
out [13. 

We have calculated the electronic structure of 
Lai.^Sr^TiOa using the LDA+DMFT method. A cu- 
bic crystal structure with 5 atoms per unit cell is uti- 
lized which is a simplified version of a fully distorted 20 
atoms/cell super lattice. Since the self-energy effects are 
crucial for the states near the Fermi energy, we treat cor- 
relations only on the downfolded t 2g orbitals of Ti atoms 
as suggested previously 0]. The Anderson impu- 
rity model is solved using Quantum Monte Carlo method 
with Hubbard parameter U = 6 eV at T = 1/(3 = 1/32 
of Ti t 2g bandwidth which delivers the self-energy E(w) 
for these orbitals using the self-consistent DMFT frame- 
work. The applicability of QMC is justified since temper- 
ature in our simulation is well below the coherence en- 
ergy, which is about 1/8 of the bandwidth. We also limit 
our consideration by dopings x larger than 10 per cent 
to stay below the coherence temperature. Once the self- 
energy is obtained, we upfold it back into the full orbital 
space which delivers the one-electron spectrum of the 
system with correlation effects taken into account. De- 
tailed description of downfolding/upfolding procedures 
to get the self-energy is given in Appendix ^ 

To treat doping away from x = the self-energy is al- 
lowed to change self-consistently while the one-electron 
Hamiltonian is assumed to be independent on doping. 
We then evaluate the frequency-dependent eigenvalues 
£kjw, ip~kju> as functions of doping. This allows us to eval- 
uate the energy and doping dependent optical conductiv- 
ity integrals both in k- and e-spaces. The integrals over 
momentum are taken on the (10 , 10, 10) mesh using the 
tetrahedron method of Ref. 131. To check the conver- 



gence we also performed the calculations on the (6,6,6) 
mesh which produces the conductivity within 5 per cent 
of accuracy. The energy integration mesh was chosen to 
have a step equal to 0.01 eV. We also broaden the imag- 
inary part of the self-energy for non-interacting bands 
with 0.0004 eV. This reproduces the LDA density of 
states of the studied compound within the accuracy of 
1-2 per cent. 

We first discuss the undoped case with x = which 
corresponds to the insulator with a small gap equal to 
0.2-0.5 eV. Model calculations for three fold degenerate 
Hubbard model, used to get the self-energy for Ti t 2g 
bands, produce a Mott-Hubbard gap equal to 2.8 eV 
but once upfolded into the LDA Hamiltonian one needs 
to take into account La 5d states in the vicinity of the 
Fermi level. The gap between the lower Hubbard band 
and La 5d bands is indeed the charge transfer gap and it 
is equal to 0.2-0.5 eV for the undoped compound. Optical 
transitions from the lower Hubbard band to La 5d give 
the main contribution to the optical conductivity in pure 
LaTi0 3 . 

Upon doping, carriers are introduced, and the system 
exhibits metallic behavior. Fig. ^ shows low-frequency 
part of <7 xx (u>) at dopings x —0.1, 0.2, and 0.3. The op- 
tical conductivity exhibits a Drude peak whose strength 
is increased with doping. The contribution to a xx {ui) 
at these frequencies is due to transitions from i) the co- 
herent part of the spectrum near the Fermi level to the 
upper Hubbard and Lanthanum bands, ii) the transitions 
from the lower Hubbard band to the upper Hubbard band 
and Lanthanum bands and iii) transitions from the lower 
Hubbard band to the coherent part of the spectra. This 
trend correctly reproduces the optical absorption exper- 
iments performed for Lai-^Sr^TiOs lSI. Comparison 
of our data with these measurements is shown in Fig. ^ 
where the measured optical conductivity at the doping 
level x = 0.1 is plotted by symbols. Overall good agree- 
ment can be found for the frequency behavior of the the- 
oretical and experimental curves. 

The strength of the Drude peak is only slightly over- 
estimated by the present theory as well as some residual 
discrepancy is seen in the region near 1 eV. We must 
emphasize that corresponding calculations based on the 
local density approximation would completely fail to re- 
produce the doping behavior due to the lack of the in- 
sulating state of the parent compound LaTi03. As a re- 
sult, the LDA predicts a very large Drude peak even for 
x — 0, which remains little changed as a function of dop- 
ing. In view of these data, the correct trend upon dop- 
ing captured by the present calculation as well as proper 
frequency behavior can be considered as a significant im- 
provement brought by this realistic DMFT study. 

More insight can be gained by comparing the effective 
number of carriers participating in the optical transitions 
which is defined by N e ff(uj c ) — J^" a(ui)dui, where m 
is free electron mass and u> c is the cut-off energy. Exper- 
imental data for N e ff{u c ) are available for the frequency 
u> c = 1.1 eV [l^. They are shown in the inset to Fig. [I] 
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FIG. 1: Low frequency behavior of the optical conductivity 
for Lai-iSr^TiOs at x = 0.1,0.2,0.3 calculated using the 
LDA+DMFT method. Experimental results |m| are shown 
by symbols for the case x — 0.1. In the inset the effective 
number of carriers is plotted as a function of doping. Squares 
show the results of the LDA+DMFT calculations. Circles 
denote the experimental data from Ref. [l9|. 



where we plot the effective number of electrons as a func- 
tion of hole concentration both from the theory and ex- 
periment |l9j . At zero doping the system is an insulator 
which gives very small N e f f for x = (this value is non- 
zero since we took uj c larger than the optical gap of the 
insulator). Upon doping, increase in N e ff is expected 
and its values as well as slope dN e f f/dx agree well with 
experiment. 

The main effect introduced by the DMFT calculation 
on the strength of the optical transitions can be un- 
derstood by looking at the Drude and interband contri- 
butions separately and comparing them with the corre- 
sponding LDA values. The LDA data give a very large 
N e ff — 1.15 which by ninety per cent consists of the 
Drude contribution. The latter can be found from the 
following equation: NRj — "^^r-f-, where plasma fre- 
quency uo p — 4.87 eV is obtained from LDA calcula- 
tions. This result is not surprising since in LDA the ti g 
states crossing the Fermi level are filled with one elec- 
tron which gives an estimation for the effective number 
of electrons participating in optical transitions at this fre- 
quency range. Thus, due to proximity to the insulator 
the DMFT suppresses ninety per cent of the Drude part 
accounted for incorrectly by the metallic LDA spectrum. 

Now we discuss optical conductivity for the frequency 
interval from to 16 eV. Fig. [21 shows a xx (io) at doping 
x = 0.1 where we compare our DMFT and LDA calcula- 
tions with the measurements in Ref. E Sharp increase 
in optical conductivity is seen at u) ~ 4 eV. This can 
be attributed to the transitions from the oxygen p-band 
into unoccupied d-states of Ti. The main peak of op- 
tical transitions is located between 5 and 10 eV which 
is predicted by both DMFT calculation (solid line) and 
the LDA (dashed line). It is compared well with the 



measured spectrum (dashed line with symbols). Since 
the self-energy corrections modify only the states near 
the Fermi level, we do not expect DMFT spectrum to 
be essentially different from the LDA one in this fre- 
quency range. Overall, the agreement at high frequen- 
cies is quite good which demonstrates reliability of the 
present method. 




Energy (Ev) 

FIG. 2: Calculated using the DMFT optical conductivity 
spectrum for La^Tii-^Os with x — 0.1 at large frequency 
interval (solid line) as compared with the experimental data 
(dashed line with symbols). The results of the LDA calcula- 
tions are shown by dashed line. 

As an additional check of the DMFT calculation, we 
have extracted the values of the linear specific heat coef- 
ficient 7 as a function of doping. Our comparisons with 
the experiment [20j are given in Fig. [3] For example, 
at x = 0.1, experimental 7 = 11 while DMFT pro- 

duces 7 equal to 14 i^rm. ■ Note that the LDA value 
here is only about 4 m " ; ^ 2 . Since DMFT renormalizes 
the density of states at the Fermi level, 7 obtained by 
this theory clearly indicates the importance of band nar- 
rowing introduced by correlations. 



IV. CONCLUSION 

In conclusion, we have shown how the optical proper- 
ties of a realistic strongly correlated system can be com- 
puted using recently developed DMFT based electronic 
structure method. We have developed a numerically 
tractable scheme which is reduced to evaluating dipole 
matrix elements as well as integrating in momentum and 
frequency spaces similar to the methods developed for 
non-interacting systems. As an application, we have 
studied the optical conductivity of Lai-^Sr^TiOa and 
found its correct dependence as a function of frequency 
and doping in comparison to the experiment. Our results 
significantly advance studies based on static mean-field 
approximations such as LDA. 

The framework that we presented should be a good 
starting point for including vertex corrections. Local ver- 
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FIG. 3: Comparison of the linear coefficient of specific heat, 7, 
as a function of doping obtained from DMFT (solid line with 
stars) and LDA (solid line) calculations against experimental 
result 20]. Experimental points are given by cross symbols 
and dot-dashed line is used as a guide for eye. 



tex corrections can be evaluated within DMFT Q while 
non-local ones can be incorporated by extending the cal- 
culations of Ref. [2lJ to the optical conductivity. This is 
analogous to how LDA spectra can be improved via the 
GW method 0. 
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APPENDIX A: COMPUTATION OF THE 
SELF-ENERGY, DMFT DOWNFOLDING AND 
UPFOLDING 

The approach described in Section ^ requires evalu- 
ation of the self-energy operator in Eq. JQ) using the 
LDA+DMFT method The latter exploits the local- 
ity of the self-energy in some orbital space, and the re- 
striction of the Coulomb interaction to a limited set of 
localized (or heavy) orbitals to be denoted by h. The rest 
of the orbitals are taken to be uncorrected (light) and 
are denoted by I. 

Notice that the locality of the self-energy is a ba- 
sis dependent statement. Under a change of the basis 
the Kohn-Sham Hamiltonian, H^, is transformed into 
UkHU^, with [/k being a unitary transformation. The 
self-energy transforms like the Hamiltonian, however, if 
£(w) is momentum independent in one basis, then in the 
new basis £' = U\J1{uj)U^ in general becomes momen- 
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turn dependent. Hence, we need to work in a very local- 
ized basis, such as the non-orthogonal LMTO's, where 
the DMFT approximation is most justified. 

Introduction of a basis set allows the partition of 
the double-counting subtracted Kohn-Sham Hamilto- 
nian H® h = Hhh — ^dc and of the Green's function into 
the light and heavy blocks: 



G(k,u) = 



n hh n hl 

ffO TT 
H lh H U / k 



Ohh Ohi 
Oi h On 

EwiM o 




(Al) 



where [•••] means matrix inversion, [i is the chemical 
potential and O is the overlap matrix. Given that the 
self-energy is local, it can be obtained from the Anderson 
impurity model 

S mp = £ 4(T)g - X a ,(T,T')c a ,(T') (A2) 
olol' ,tt' 

+ E ^c+(r)4(r)c 7 (T) Cj (r), 

af3~fS,T 

where Qq is the bath Green's function which obeys 
the self-consistency condition [B| generalized to non- 
orthogonal basis set: 



(w + /x)O-iJ (k) -E(w) 



hh 



(A3) 

When a group of bands is well separated from the oth- 
ers it is possible to recast the previous self-consistency 
condition at low frequencies in a form which resembles 
the DMFT equations derived from a Hamiltonian involv- 
ing the h degrees of freedom only. In the one-electron 
approach it goes under the name downfolding p^ . 

Performing standard matrix manipulations and a low- 
frequency expansion with linear accuracy in u> (which is 
justified for low-energy calculations provided the separa- 
tion of energy scales between the band near the Fermi 
level and the rest) we rewrite the heavy block of the 
Green's function as: 



G hh (k,uj)= Z^u - H{k) - £ 



h h 



(A4) 



where renormalization amplitude and effective Hamil- 
tonian are given by 



z: 



O hh +K hl K^O u K^K lh 
-OuK^Kin - KuK^Oih, 



H(k) 



HI 



hh 



ih- 



(A5) 



Here 7 stands for a pair of indices / or h. Finally we 
perform a unitary transformation S in the heavy block, 



so as to work in a nearly orthogonal basis in the /i-sector: 



s^[J2z k f 1 s = i. 



(A6) 



Applying this transformation to Eq. I|A5J| we arrive to 
the local Green's function in the new basis 

G hh (co) =£> + M)O e //(k) - H eff (k) - EHP 1 , 

(A7) 

and to a new DMFT self-consistency condition: 

0mh(") = Gm ( A8 ) 

This set of equations has clearly the form of the DMFT 
equations of a model involving heavy electrons only, with 
a Hamiltonian and an overlap matrix: 

O eff (k) = S^Z^S, (A9) 
H eff (k) = S^H(k)S + f iO eff (k). (A10) 

The self-energy E is still computed from the Ander- 
son impurity model, but the Coulomb interaction of this 
model is renormalized to a smaller effective interaction 
U e f f matrix 



U 'eff,aif3>y5> 



(All) 

Until now the discussion is general, and applies to any 
system where there is a set of bands well separated from 
the rest. Further simplifications are possible, if we as- 
sume that the system has cubic symmetry and that the 
overlap O e // is the unit matrix. For d-electrons, cu- 
bic symmetry makes the self-energy and local Green's 
function diagonal. In this case the momentum sum in 
Eq. i|A7|l can be replaced by the integral over energy. 
The local Green's function can be calculated as a Hilbcrt 
transformation 



+ OO 



de ■ 



D(e) 



+ PL — 



(A12) 



Here, D(e) is the density of states of the the reduced 
Hamiltonian H e ff(k). Notice that the cubic symmetry 
keeps U e ff diagonal if the bare Coulomb matrix U has 
that property. 

Upfolding is a procedure which is "inverse" to the 
downfolding described above. One simply converts the 
self-energy S obtained from the DMFT calculation into 
the block self-energy = ST,S> , which is to be in- 
serted to the original LDA Hamiltonian, in order to com- 
pute the local Green's function G(u>). 

In general, the downfolded density of states D(e) ob- 
tained from H e f f has a non-zero first energy moment and 
depends in a non-linear way on the value of the double 
counting correction, as well as on the chemical poten- 
tial which enters the formulation of the original problem 
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containing all electronic bands. Furthermore, the value of 
the chemical potential in the LDA+DMFT calculations 
does not need to be the same as the LDA value. 

The reduction of the self-consiste nt LD A+DMFT 
equations to the form described by Eq. (|A12I) with D(e) 
being the partial LDA density of states of the heavy or- 
bitals was suggested and used in Ref. • Unfortunately, 
this partial density of states contains weight at high ener- 
gies, and if this is omitted, the normalization condition is 
violated. The derivation presented in this Appendix elim- 
inates these difficulties, and instead suggests an alterna- 
tive procedure in which we first carry out a tight -binding 
fit of the LDA bands (downfolding) near the Fermi level, 
and then use it to estimate D{e). Our derivation also in- 



dicates how one goes back (i.e. upfolds the self-energy) 
to the all-orbital Hamiltonian. In our calculations us- 
ing the downfolded equations fi was adjusted to get the 
correct density of d-electrons. In the upfolded Green's 
function fj, was taken to be the LDA chemical potential, 
and Srfc was deduced from a constant shift of the heavy 
orbitals by obtaining the total number of electrons from 
the integral of the spectral function 

k a/3 

multiplied by the Fermi function. 



